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ABSTRACT 



Context. When type la supernovae (SNe la) were chosen as distance indicator to measure cosmological parameters, Phillips relation 
was applied. However, the origin of the scatter of the maximum luminosity of SNe la (or the variation of the production of 56 Ni) is 
still unclear. The metallicity and the carbon abundance of white dwarf (WD) before supernova explosion are possible key, but neither 
of them has an ability to interpret the scatter of the maximum luminosity of SNe la. 

Aims. In this paper, we want to check whether or not the carbon abundance can be affected by initial metallicity. 
Methods. We calculated a series of stellar evolution. 

Results. We found that when Z < 0.02, the carbon abundance is almost independent of metallicity if it is plotted against the initial 
WD mass. However, when Z > 0.02, the carbon abundance is not only a function of the initial WD mass, but also metallicity, i.e. for 
a given initial WD mass, the higher the metallicity, the lower the carbon abundance. Based on some previous studies, i.e. both a high 
metallicity and a low carbon abundance lead to a lower production of 56 Ni formed during SN la explosion, the effects of the carbon 
abundance and the metallicity on the amount of 56 Ni are enhanced by each other, which may account for the variation of maximum 
luminosity of SNe la, at least qualitatively. 

Conclusions. Considering that the central density of WD before supernova explosion may also play a role on the production of 56 Ni 
and the carbon abundance, the metallicity and the central density are all determined by the initial parameters of progenitor system, i.e. 
the initial WD mass, metallicity, orbital period and secondary mass, the amount of 56 Ni might be a function of the initial parameters. 
Then, our results might construct a bridge linking the progenitor model and the explosion model of SNe la. 

Key words. Stars: white dwarfs - stars: supernova: general 



1. Introduction 

In their function as one of the distance indicators, type la super- 
novae (SNe la) show their importance in determining cosmolog- 
ical parameters, which resulted in the discovery of the acceler- 
ating expansion of the universe (Riess et al. 1998; Perlmutter et 
al. 1 1 999l >. The result was exciting and suggested the presence of 
dark energy. At present, SNe la are proposed to be cosmological 
probes for testing the evolution of the dark energy equation of 
state with time and testing the evolutionary history of the uni- 
verse (Riess et al. 2007, Kuznetsova et al. 120081 Howell et al. 
I2009ab . 

When SNe la are applied as a distance indicator, the Phillips 
relation is adopted, which is a linear relation between the ab- 
solute magnitude of SNe la at maximum light and the magni- 
tude drop of the B light curve during the first 15 days follow- 
ing the maximum (Phillips 1993). This relation was motivated 
by the observations of two peculiar events, i.e. SN 1991bg and 
SN 199 IT, and implies that the brightness of SNe la is mainly 
determined by one parameter. There is a consensus that a SN 
la is from a thermonuclear explosion of a carbon-oxygen white 
dwarf (CO WD) and the amount of 56 Ni formed during the su- 
pernova explosion dominates the maximum luminosity of SNe 
la (Arnett 1982). But the origin of the variation of the amount 
of 56 Ni for different SNe la is still unclear (Podsiadlowski et 
al. 12008). Many efforts have been paid to resolve this problem. 
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Some multi-dimensional numerical simulations showed that the 
ignition intensity (the number of ignition points) in the cen- 
ter of WDs or the transition density from deflagration to deto- 
nation are wonderful parameters interpreting the Phillips rela- 
tion (Hillebrandt & Niemeyer |2000l Hoffich et al. 120061 120101 
Kasen et al. 120101) . In addition, the ratio of nuclear-statistical- 
equilibrium (NES) to intermediate-mass elements (IME) in the 
explosion ejecta is likely the key parameter for the width of SN 
la light curve and its peak luminosity (Pinto & Eastman 2001 ; 
Mazzali et al. 200Q 120071) . For the simulations above, a rea- 
sonable assumption is that these parameters are determined by 
one or some properties of SNe la progenitor. Since these pa- 
rameters are free ones for the numerical simulations of SNe la 
explosion, the question was turned into which property or prop- 
erties of progenitor system determine these parameters. This is 
still unclear. Lesaffre et al. (2006) carried out a systematic study 
of the sensitivity of ignition conditions for H-rich Chandra sin- 
gle degenerate exploders on various properties of the progeni- 
tors, and suggested that the central density of the WD at ignition 
may be the origin of the Phillips relation (see also Podsiadlowski 
et al. 120081) . This suggestion was uphold by detailed multi- 
dimensional numerical simulations of explosion (Rrueger et al. 
2010). But in the models of Hoflich et al. (2010), the central 
density is only a second parameter, and the cooling time of the 
WDs before mass transfer in Lesaffre et al. (2006) an d Rrueger 
et al. ( 120101) is shorter than 1 Gyr. However, there are SNe la as 
old as 10 Gyr. The WDs with such a long cooling time may be- 
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come more degenerate before the onset of accretion phase. Some 
other processes like C and O separation or crystallization may 
occur, and dominate the properties of the CO WD (Fontaine et al. 
1200 11 1. How the extremely degenerate conditions affect the prop- 
erties of SNe la is still unclear. Then, the suggestion of Lesaffre 
et al. ( 2006 ) should be checked carefully under extremely degen- 
erate conditions (Bravo et al. 2010b >. Some numerical and syn- 
thetical studies showed that metallicity has an effect on the final 
amount of 56 Ni, and thus the maximum luminosity (Timmes et 
al. 12003] Trav aglio et al. 120051 Podsiadlowski et al. l2TjQ6l Bravo 
et al. 120 1 0al > and there do be some observational evidence of 
the correlation between the properties of SNe la and metallicity 
(Branch & Berg h[T993l Hamuy et al. [19961 Wang et al. [19971 
Cappellaro et al l 19971 Shanks et al. l2002l l. However, the metal- 
licity seems not to have the ability to interpret the whole scat- 
ter of the maximum luminosity of SNe la (Timmes et al. 120031 
Gallagher et al. [20081 Howell e t al. |2009bl l. 

Nomoto et al. (119991 120031 1 suggested that the average^] ra- 
tio of carbon to oxygen (C/O) of a white dwarf at the moment 
of explosion is the dominant parameter for the Phillips relation 
(see also Umeda et al. [1999b ). The higher the C/O, the larger the 
amount of nickel-56, and then the higher the maximum luminos- 
ity of SNe la. By comparing theory and observations, the results 
of Meng, Chen & Han ( f20D9b a nd Meng & Y ang d2010al) upheld 
this suggestion. Nomoto et al. dl999l [2003 > used total 12 C mass 
fraction included in the convective core of mass M = 1.14M 
just before the SN la explosion, X(c), to represent the C/O ratio 
and their suggestion is based on a fact that the dependence of 
X(c) on metallicity is small when it is plotted against the initial 
mass of WDs when Z < 0.02. Actually, when Z < 0.02, the ef- 
fect of metallicity on the amount of 56 Ni can be neglected. Only 
when Z > 0.02, metallicity has a significant influence on the pro- 
duction of 56 Ni in SNe la explosion (Timmes et al. I20031 >. Here, 
we want to check whether or not the C/O ratio still can not be 
affected by metallicity when Z > 0.02. 

In section[2l we simply describe our method, and present the 
calculation results in section[3] In section|4] we show discussions 
and our main conclusions. 

2. METHOD 

In this paper, our work is based on the single degenerate sce- 
nario, i.e. a SN la is from a CO WD in a binary system and 
its companion is a main-sequence or a slightly evolved star 
(WD+MS), a red giant st ar (WD +RG) or a helium star (WD 
+ He star) (Whelan & Iben l 19731 Nomoto, Thielemann & Yokoi 
1984). This scenario has been widely studied by many groups 
(Li & van den Heuvel [T997l Hachisu et al. U999al Langer et al. 
[20 001 Ha n & Podsiadlowski| 2004l C hen & L il2007l Hachisu et 
al. 120081 Men g, Chen & Han |2009l Lii et al.|2009l Wang et al. 
I2009al I2009bl Meng & Yang 1201 Obi 1201 Oct Wang, Li & Han 
2010). We assume that an initial CO WD is derived from a main 
sequence (MS) star in a primordial binary system (primary). The 
CO WD accretes hydrogen-rich material from its companion via 
Roche lobe overflow or wind, where the companion is a normal 
star. The accreted hydrogen-rich material is burned into helium, 
and then the helium is converted to carbon and oxygen. The CO 
WD increases its mass until the mass reaches 1.378 M (close 
to the Chandrasekhar mass limit, Nomoto, Thielemann & Yokoi 
1984) where it explodes in a thermonuclear supernova. A binary 
systems with the same primordial primary but different orbital 

1 If there is no special statement, the C/O ratio in this paper is the av- 
erage value over the whole WD structure just before the SN explosion. 



period may produce a CO WD with different mass. For simplic- 
ity, we assume that if the MS mass of the primary in a primordial 
binary system is same, the initial mass of the CO WD from the 
binary system is same and is equal to the core mass derived from 
the envelope-ejection model in Han, Podsiadlowski & Eggleton 
( 1 19941 ) and Meng et al. (2008 ) (see below in details). Because the 
primary in a binary system may lose its hydrogen envelope due 
to the influence of secondary before it fulfills the criterion for 
the envelope ejection, the initial mass of CO WD in this paper 
should be taken as an upper limit for a real case. Actually, this 
assumption can not affect our results since we only want to find 
a relation between the C/O ratio, the WD mass and metallicity. 
The method used here is similar to that in Umeda et al. (|1999a) 
and Hoflich et al. d20T0b . 

The C/O ratio before SN la explosion is a result of stellar and 
binary evolution, i.e. during central helium burning and thin shell 
burning during the stellar evolution and the accretion to close to 
the Chandrasekhar mass limit (1.378 M , Nomoto, Thielemann 
& Yokoi 1984). After the central helium burning phase, the C/O 
ratio is low, i.e. 0.25 - 0.5 depending on the initial mass and 
metallicity of main sequence (MS) star (Umeda et al. [1999a). 
The C/O ratio obtained from the burning shell is » 1, because the 
helium in the shell has a lower density and a higher temperature 
compared to helium burning in the core (Hoflich et al. 12010) . 

In the paper, we calculate a series of stellar evolutions with 
the primordial MS mass from 1 M G to 6.5 M until the stars 
evolve to the asymptotic giant branch (AGB) stage. When a star 
evolves to the stage, its envelope may be lost if the binding en- 
ergy (BE) of the envelope transforms from a negative phase to 
positive one (Paczyriski & Ziolkowski l 19681 see also Fig. 2 in 
Meng et al. 120081 ). We calculated the BE of the envelope by 

f M * Gm 

AW= ( + U)dm, (1) 

Jm c r 

where M c is the core mass, M s is the surface value of the mass 
coordinate m. U is the internal energy of thermodynamics where 
those due to ionization of H and dissociation of H2, as well as 
the basic |2i77ju for a perfect gas are all included). Here, we 
assume that a star will lose its envelope when the BE of the star's 
envelope increases to the point of AW = and the core mass 
at the point is the final WD mass. The method here is robust 
and its virtue is significant because we need not consider the 
specific mechanism of mass loss since the mass loss rate is very 
uncertainty (see Meng et al. 2008 in details about this method). 
We assume that the remnant after envelope ejection is a CO WD 
if carbon and oxygen have not been ignited at the moment of 
envelope ejection. Following shell burning, we assume that the 
C/O ratio is 10 until M WD = 1 .378M as did in Umeda et al. 

dl999al ) and Hoflich et al. d20T0b . 

We use the stellar evolution code of Eggleton (119711 119721 
119731) . which has been updated with the latest input physics over 
the last three decades (Han, Podsiadlowski & Eggleton 1994; 
Pols et al. 19951 1 19981 ). The chemistry of a WD is mainly deter- 
mined by the competition of two major nuclear reactions pow- 
ering the He burning, i.e. 3a and l2 C(a, y) le O. As discussed by 

2 In fact, for the core chemical profile produced during the AGB 
phase, the steady increase of the carbon abundance in the region ad- 
jacent to the inner flat profile is synthesized by the shell in the early 
AGB, whereas the following abrupt peak is left by the thermally pulsing 
phase. The slop of the profile produced during the early AGB depends 
on both the initial mass and metallicity, and the C/0= 1 is not always 
a good approximation. However, to compare with previous results, we 
still simply make a constant-ratio assumption. 
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Imbriani et al. (1200 11 1 and Prada Moroni & Straniero (2002), the 
final C/O after the central helium exhaustion not only depends on 
the rate of these two reactions, but also is significantly influenced 
by the efficiency of convective mixing operating at the central 
helium burning phase. We set the ratio of mixing length to local 
pressure scale height, a = l/H p , to 2.0, and set the convective 
overshooting parameter, <5ov, to 0.12 (Pols et al. 1997, Schroder 
et al. 1 19971 ), which roughly corresponds to an overshooting 
length of 0.25Hp. The two parameters are adopted during the 
whole evolution of star. For the occurrence of convective insta- 
bility near the He exhaustion in the central core, some breathing 
pulses are expected. (Castellani 1985 1989). This phenomenon 
may also occur naturally in our code. The reaction rates are from 
Caughlan & Fowler ( 1988 ), except for the l2 C(a, y) 16 reaction 
which is taken from Caughlan et al. (1985 ). The range of metal- 
licity is from 0.0001 to 0.1, i.e. 0.0001, 0.0003, 0.001, 0.004, 
0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.08, 0.1. The o pacity tables 
for these metallicties are compiled by Chen & Tout (2007 1 from 
Iglesias & Rogers ( 1996 ) and Alexander & Ferguson ( 1994 1. For 
a given Z, the initial hydrogen mass fraction is assumed by 



X = 0.76 - 3.0Z, 



(2) 



(Pols et al. 119981 ), and then the helium mass fraction is Y = 1 - 
X - Z = 0.24 + 2Z. Based on the correlation between X, Y and 
Z used here, Pols et al. ( 1998 ) accurately reproduced the color- 
magnitude diagrams (CMD) of some clusters. 

In this paper, we skip thermal pulses by taking a longer time- 
step to reduce computing time, i.e we study the average evolu- 
tion of thermally pulsing AGB models. This treatment on the 
thermal pulses may lose some information about the structural 
change of the envelope due to thermal pulse. However, these 
treatment will not affect our final results because the situation 
of AW = is fulfilled before the onset of thermal pulse and after 
the point a simple assumption of C/O — 1 is adopted (see also 
in Meng et al.|2008). 

3. Results 

In this paper, we use total carbon abundance to represent the C/O 
ratio and a core mass to represent the initial WD mass as did in 
Umeda et al. (1999b). In figure [1] we show the abundances of 
several elements in mass fraction in the inner core of three rep- 
resentative models with different metallicity at the moment of 
envelope ejection. In the figure, the H burning shell is located 
around the mass coordinate of m = O.9M , which means that 
the remnants from the models have a similar mass, i.e. the ini- 
tial WD mass from the stars is similar. An interesting feature in 
the figure is that the carbon abundance and the flat part of the 
carbon abundance profile in the inner core decrease with metal- 
licities (Umeda et al. [T999a; Dominguez et al. 2001). The car- 
bon abundance is mainly derived from the result of the competi- 
tion between 3a and l2 C(a, y) ie O reaction in the central helium 
burning phase, while the flat profile is the result of the central 
helium burning, which occurs in a convective core. An increase 
of the metallicity leads to an increase of the radiative opacity, 
and consequently, a decrease of the central temperature at given 
phase, the decrease of the He core mass at the 3a onset. For the 
central helium burning phase, this leads to a smaller convective 
core, hence a smaller region characterized by a flat C/O profile. 
A low central temperature favors the destruction reaction of I2 C, 
namely the l2 C(a, y) I6 with respect the 3a reaction, which is 
the main cause for the decrease of the central carbon abundance. 
In addition, the central carbon abundance is also relevant to the 
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Fig. 1. Abundances of several elements in mass fraction in the 
inner core of three models with different primordial MS masses 
and different metallicities at the moment of envelope ejection. 
The solid, dotted, dashed and dot-dashed lines represent hydro- 
gen, helium, carbon and oxygen abundance, respectively. Top: 
Z = 0.001 and Mi = 3.51; Middle: Z = 0.02 and Mi = 4.52; 
Bottom: Z = 0.06 and M { = 4.32. 



helium abundance in central helium burning phase, and then to 
metallicity via Y — 1 - Z at the beginning of central helium 
burning. For a given temperature, a low helium abundance, i.e. 
a high metallicity, means a slight lower burning rates of 3a re- 
action. However, between the model of Z = 0.001 and Z = 0.02 
in figureQ] the difference of carbon abundance is not significant, 
but the difference between the model of Z = 0.02 and Z = 0.06 is 
remarkable. Then, the effect of metallicity on the C/O ratio may 
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Fig. 2. The relation between the total carbon abundance of WDs 
before supernova explosion and the initial WD mass for SNe la 
with different metallicities. The gray points represent those that 
may not contribute to SNe la. 
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Fig. 3. The relation between the total carbon abundance of WDs 
before supernova explosion, the initial WD mass and metallicity. 
The solid line is the best fitted linear relation, where the line 
fits all the points shown in Fig. [2] The lower panel shows the 
difference between the points and the fitted line. 

not be neglected when Z > 0.02. In addition, the hydrogen abun- 
dance decreases and helium abundance increase with metallicity 
in the figure, which is a natural result of equation |2fl 

In figure [Jfl, we show the relation between the total carbon 
mass fraction of WDs before SNe la explosion and the initial 
WD mass for SNe la with different metallicities. For the cases 
of Z < 0.02, the results here is similar to that in Nomoto et 
al. ([T999, 2003 ), i.e. the relation between the carbon abundance 
and the initial WD mass is independent of metallicity, especially 
for M WD > O.8M . When M WD < O.8M , the scatter is en- 
larged for the cases of Z < 0.02, which is mainly from low- 

3 The hydrogen abundance and the helium abundance are not accu- 
rately equal to those derived from equation [2] because of the first and 
second dredge-up (Busso et al. H999t . 

4 Here, the gray points are plotted based on the study in Meng et al 
(2009 ). Please keep in mind that the results for the lower limit of WDs 
for SNe la in Meng et al. ( 2009 1 may be higher than that in Meng & 
Yang (2010a) since some special effects such as mass-stripping effect 
of an optically thick wind and the effect of thermal unstable disk are 
considered in Meng & Yang d2010ab . 
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Fig. 4. The relation between the total carbon abundance of WDs 
before SNe la explosion, the initial WD mass and metallicity, 
where the points which may not contribute to SNe la are cut off. 
The solid line is the best fitted linear relation. The lower panel 
shows the difference between the points and the fitted line. 



metallicity cases (< 0.001, see also Fig. 12 in Umeda et al. 
1999a). However, for the low-metallicity cases, the WDs with 
a mass smaller than O.8M may not contribute to SNe la (see 
figure 5 in Meng, Chen & Han|2009|. 

A remarkable feature in figure |2]is that when Z > 0.02, the 
relation between the carbon abundance and the initial WD mass 
significantly deviates from those of Z < 0.02 and the deviation 
increases with metallicity, i.e. for a given initial WD mass, a high 
metallicity leads to a lower carbon abundance as shown in figure 
Q] This is similar to that found by Timmes et al. (2003 ), i.e. only 
when Z > 0.02, the influence of metallicity becomes significant. 
So, the discovery of Umeda et al. (1999b) and Nomoto et al. 
( 1999, 2003 ) is a low-metallicity limit of the result found in this 
paper. 

In the three-dimensional space of (Xc, M core , Z), the relations 
between the carbon abundance and the initial WD mass with 
different metallicity are almost in a plane. We take the Xc as 
a function of M core + aZ and use the minimum x 2 method to 
find the best relation. We found that a linear relation may well 
represent the relation when a = 3.5. The best fitted relation is 

X c = 0.5531 - 0.2528(M core + aZ), (3) 

where = 5.16 x 10~ 2 . We also try to use a parabola to fit 
the relation between Xc and M cole + aZ, but the improvement is 
slight, i.e. x^tAn ~ 4-84 x 10" 2 . Seen from the equation [3] both 
a high initial WD mass and a high initial metallicity lead to a 
lower carbon abundance. In addition, if the metallicity is low 
enough, its effect on the carbon abundance can be neglected. But 
the effect of a high metallicity is significant. For example, for a 
typical value of M\ = O.8M (Meng, Chen & Han 120091 Meng 
& Yang 2010a ), the uncertainty of Xc derived from metallicity 
for Z < 0.02 is less than 5%, but may be as large as 25% for 
Z > 0.02. 

If the points which may not contribute to SNe la in Fig. [2] 
are cut off, we even may obtain a more tight linear relation, i.e. 
X c = 0.5824 - 0.2862(M core + aZ) where a = 3.0 and x 2 min = 
1.46 x 10" 2 (see Fig. gj). 
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4. Discussions and conclusions 

When one simulates the SN la explosion, the carbon abundance 
(the C/O ratio) and metallicity are always taken as free parame- 
ters. But neither of them has an ability to interpret the variation in 
the mass of 56 Ni ejected by SNe la, especially for subluminous 
SNe la (such as 1991bg-like supernovae, Timmes et al. 2003; 
Ropke et al. 12006a; Bravo et al. I2010bl l. However, many evi- 
dence shows the correlation between the maximum luminosity 
of SNe la and metallicity. For example, many groups noticed 
that subluminous SNe la occur exclusive in massive galaxies 
(Neill et al. 120091 Su llivan et al. 120101 Gonzalez-Gaitan et al. 
120101 Lampeitl et al. 1201 01) . Considering the mass-metallicity 
relation of galaxies (Tremonti et al. 12004) . subluminous SNe la 
favor more metal-rich environments. In addition, some 1991T- 
like SNe la were discovered in metal-poor environments^ (Prieto 
et al. 120081 Badenes et al. 120091 Khan et al. 120101) . So, metallic- 
ity should play a more significant role in a rather wide range than 
that suggested in theory. The connection between theory and ob- 
servations should be resolved. 

In this paper, we found that when Z > 0.02, the carbon 
abundance before SN la explosion is affected by not only ini- 
tial WD mass but also metallicity. For a given initial WD mass, 
the higher the metallicity, the lower the carbon abundance. The 
relation between the carbon abundance, the initial WD mass and 
the metallicity may be represented by a simple linear relation. 
Thus, when one simulates the SN la explosion, the carbon abun- 
dance (the C/O ratio) and metallicity would not be taken as free 
parameters. Their effect on the amount of 56 Ni is enhanced by 
each other since the effect of a high metallicity and a low car- 
bon abundance on the amount of 56 Ni is similar, i.e. producing a 
lower amount of 56 Ni (Nomoto et al.[[999 2003; Timmes et al. 
2003). So, our result may account for the variation of the max- 
imum luminosity of SNe la, at least qualitatively providing a 
method to conquer the confliction stated above. However, please 
keep in mind that whether or not the combined action of the C/O 
ratio and metallicity has enough ability to interpret the variation 
of the maximum luminosity of SNe la still should be verified 
carefully. So, we encourage someone to do a detailed numerical 
simulation on this problem based on the result here. 

Furthermore, the central density at the moment of supernova 
explosion may also play a role to some extent (Hoflich et al. 
2010; Krueger et al 1201 Oi l, and the carbon abundance, metal- 
licity and the central density all contribute to the variation of 
the maximum luminosity of SNe la (Ropke et al. 1200 6b). By a 
simple assumption that the carbon abundance is the function of 
initial WD mass and the central density is determined by initial 
WD mass and its cooling time, Meng et al. (120101 1 noticed that 
the WDs with a high carbon abundance usually have a lower 
central density at ignition, while those having the highest cen- 
tral density at ignition generally have a lower carbon abundance. 
Interestingly, the effect of a high metallicity, a low C/O ratio and 
a high central density on the amount of 56 Ni is similar, i.e. pro- 
ducing a less amount of 56 Ni and thus a dimmer event (Nomoto 
et al. [19991 120031 Timmes et al. l2003l Krueger et al l20"T0l l. then 
all of them could contribute the fact that elliptical galaxies fa- 
vor the dim SNe la (Hamuy 1996). Although it is still unclear 
how the metallicity affects the central density, we may hypoth- 
esize optimistically that it is unreasonable that take the carbon 

s Even 2003fg-like SNe la which are over-luminous events favor 
metal-poor environments (Taubenberger et al. 2010 1, but the discussion 
about these SNe la is beyond the scope of this paper since these SNe la 
are over Chandrasekhar mass limit (Howell et al. 2006 1 and this paper 
is based on the Chandrasekhar mass model. 



abundance, metallicity and the central density as free parame- 
ters when one simulates SN la explosion, and these parameters 
may all contribute the production of 56 Ni (Ropke et al. [2006b), 
which could be a key to open the origin of the Phillips relation. 
Moreover, these parameters can be determined by the initial pa- 
rameters of a progenitor system, i.e. the WD mass, metallicity, 
orbital period and secondary mass. For example, the C/O ratio is 
a function of the initial WD mass and metallicity, while the cen- 
tral density of a WD before supernova explosion is determined 
by accretion rate and its cooling time before the onset of mass 
transfer, which is related with the initial WD mass, secondary 
mass and period. Then our study might provide a bridge linking 
the progenitor model and explosion model of SNe la. 
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Appendix A: Some physical quantities for models 

In this paper, we calculated a large and fine grid models. Here, 
we provide some physical quantities of the models such as CO 
core mass, total carbon abundance X(C), initial main-sequence 
mass and metallicity, which may be helpful for completing the 
resolution in the field. Here, the boundary of the CO core was 
located at the 10% of the maximum helium abundance in the 
helium shell (see the left dotted line in Fig. [T}. 
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Table A.l. The CO core masses (in M ) for different metallicities (in Z , Row 1) and different initial masses (Column 1, in M ). 
The bars represent models whose final fates are ONeMg white dwarfs. 
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Table A.2. The total carbon mass fraction of WDs before SNe la explosion for different metallicities (in Z , Row 1) and different 
initial masses (Column 1, in M ). The bars represent models whose final fates are ONeMg white dwarfs. 
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